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A model describing the coincidental timing of scintillating fibers is developed. Fiber geometry, the rate of scintillation decay 
together with the mean number, spatial dispersion and attenuation of emitted photons is considered. For a specific selection of 
probability distributions and parameters involved, the entire coincidental timing distributions, corresponding FWHM values and 
the photon detection efficiencies are extracted. The significance of the number of photons from the scintillation process is specially 
emphasized. Additionally, the model is extended to include a triggering feature, experimentally realized by coupling fibers to 
any photon resolving device. Finally, the measurements of a coincidental timing distribution were performed, with an excellent 
agreement found between the experimental and predicted theoretical results. 
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1. Introduction 

Scintillating fibers have been extensively investigated from 
both theoretical and experimental point of view. Consequently, 
their scintillation mechanisms and light transmission properties 
are well known and conveniently summarized in [1], together 
with a list of the most relevant references on the subject. Nowa- 
days, there is an active effort to include fibers as an integral part 
of sophisticated detector systems for nuclear experiments. Be- 
ing intended for a detection of charged particles and/or elec- 
tromagnetic radiation, fibers are required to provide an ade- 
quate, if not excellent coincidental timing resolution, vital for 
precise reconstruction of the scintillating pulse position along 
their length. Compared to the resolution of bulk scintillators - 
an order of magnitude below the nanosecond scale [2] - early 
measurements indicate that attainable values are well above this 
range [3]. The rate of scintillation process, together with a low 
number of detected photons being subject to a wide spatial dis- 
persion, is considered to be the main cause for the coincidental 
timing discrepancy with that of bulk scintillators. Commonly, 
Monte Carlo simulations present a natural tool for estimation 
and prediction of experimental results, easily capable of in- 
cluding all the required physical principles. However, in this 
paper an analytical model is developed, describing an asymp- 
totic form of otherwise measured coincidental timing distri- 
butions, giving rise to the central resolution defining FWHM 
value. Compared to the earlier, long-established models [4, 5], 
the one developed herein extends beyond the sole scintillation 
process, taking into account a subsequent light propagation, 
while proposing the simple manner in which to include even 
the effects of the photon detecting units. An assumption of a 
fully general form for the emission of the scintillation photons. 
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their spatial distribution and subsequent attenuation inside the 
fiber material allows physical considerations of varying com- 
plexity to be adopted - from a simple meridional approximation 
to a considerably more refined description of optical processes 
involved. Due to the level of generality acchieved, though de- 
veloped with scintillating fibers in mind, the model is indis- 
criminately applicable to any kind of light guides, making it 
an acute mathematical tool for predicting the experimental re- 
sults far beyond the assumptive limitations of a specific setup. 
As a starting point a model for the photon propagation times 
is considered - initially developed in [2] for bulk scintillators 
and later applied to scintillating fibers in [6] - from which the 
technical formalism was adopted. 



2. Basic model 

In common nuclear or particle physics experiment the arrival 
time of a signal is determined by the leading edge or the con- 
stant fraction discrimination of the signal's leading edge. Since 
the first photon impinging on the detector initiates the rise of 
the signal, its statistics is of the utmost importance for the de- 
scription of signal timing properties. 

Therefore, let us assume that inside the scintillating fiber of 
a length L at a distance / from one of its ends a scintillation 
pulse was induced (Fig. 1), emitting a total of n photons. To 
determine a probability /„ for the first arriving photon to reach 
the fiber end, several separate cases must be considered. For 
example, first arriving photon may correspond to the first emit- 
ted (probability pi). On the other hand, first m-l emitted ones 
may be lost either by escaping from the fiber before reaching 
its end or by absorption inside the fiber material. Therefore, 
m-th emitted one becomes first to reach the fiber end (probabil- 
ity /?,„) and induce the signal in the detector. Consequently, all 
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The second useful term A: 





Figure 1 : A scintillating pulse induced by the chaiged pailicle q passing through 
the fiber at a distance / from one of its ends. The difference of initial photons' 
times of arrival at each end constitutes a basis for measuring the coincidental 
timing resolution. A wide spatial dispersion of isotropically emitted photons is 
also illustrated, increasing the emission probability for greater angles, while 
degrading attainable resolution. 



such contributions give rise to probability /„: 



fn — ^ ^ Pn 
m=l 



(2.1) 



Let us beforehand define the moment for the first photon 
arrival at the fiber end, the moment of its emission and the 
time tp needed for its propagation through the fiber The fol- 
lowing holds: 

fA = f£ + tp (2.2) 

Most simplistic in form, previous relation will be fundamental 
for final calculations. For a detailed analysis, let us separate 
contributions to pm'- 



(2.3) 



with a„, regulating m - 1 photons emitted before the m-th one, 
/3m regulating m-th emitted one as first arriving to the fiber end, 
and regulating all the n-m following ones. For a description 
of these probability contributions let us assume a completely 
general form £(?£) for the emission probability distribution and 
S (tp) for the spatial distribution of emitted photons expressed 
in terms of the photon time propagation, i.e. translated into 
the path length dispersion. Furthermore, an attenuation factor 
A(tp) governing the probability of single photon reaching the 
fiber end will be required. Defining two additional terms will 
prove to be most convenient for fuither calculations. Therefore, 
for the fiber excitation occurring at the moment ?£ = 0, let us 
define !(?£): 



I(f£) 



(2.4) 



'max 

As J Sit„)Ait^)dt„ 



(2.5) 



denotes the probability for emitted photon to actually reach the 
fiber end. With f„„„ as minimal time required for the photon 
propagation' and t„,ax as maximal time of propagation permit- 
ted^, A is governed by the probability for photon to become 
trapped by internal reflections inside the fiber and not to be ab- 
sorbed within the fiber material. 

For the m-th emitted photon to become the first to aiiive at the 
fiber end, those m - 1 previously emitted have to be lost, either 
by escaping the fiber or by absorption. Therefore, considering 
the photon combinations, factor from (2.3): 



n 

m — 1 



[I(f£)(l - A)] 



m-1 



(2.6) 



is given by the probability for initial m - 1 photons to be emit- 
ted prior to the emitting moment of the m-th one, and subse- 
quently not to reach the fiber end. It is to be noted that a,,, was 
constructed without arranging the lost photons in time, which 
is an approach validated in Appendix A. Furthermore, after 
being emitted at f^, factor/?,,,: 



/3„, = (n-m+l)E{tE)A{tp)S(tp) 



(2.7) 



regulates the spatial direction, i.e. path length, and attenuation 
probability for the m-th photon. The number of remaining sin- 
gle photon selections is also taken into account. Finally, factor 

7m- 



rm = [l-I(f£)r 



(2.8) 



as a probability for photon emission prior to the moment ?£. 



is restrained only by the requirement for the emission of re- 
maining n - m photons occurring after fg, regardless of their 
outcome. Isolated combinatory factor is absent because all the 
selection options were depleted by a,,, and . It is to be noted 
that am and y,,, are true probabilities, while /3,„ is, in fact, a 
probability density. 

With am, Pm, 7m obtained, p,,, is completely determined by 
(2.3), and /„, consequently, by (2.1). Therefore, writing /„ ex- 
plicitly: 

/„ =E(tE)A(tp)S(tp)x 

X J(„ - ,„ + 1)( \ [l(tE)(l - A)]'"-' [1 - KtE)]"-"' 

m=\ ^ ^ 

(2.9) 

it may be noted that by shifting a summation index a step back- 

' Minimal propagation time f,,,,,, corresponds to the shortest distance path 
between the point of scintillation inside the fiber and the fiber end, i. e. to the 
photon emitted along the fiber axis. 

^Maximal propagation time t^ax is commonly considered to be defined by 
the optical condition for a total reflection off the fiber walls, i.e. by the criti- 
cal angle for total reflection. However, this is not strictly true, which will be 
discussed in Appendix C. 
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wards, a binomial expansion remains, reducing (2.9) into: 



3. Fine corrections 



Utp,tE\T)^nE(tE)A{tp)S,{tp)[\-\{tE)Kir-' (2.10) 

In (2.10) an explicit dependency on the position I of a scintil- 
lation pulse origin along the fiber, the emission moment f£ and 
the photon propagation time tp was written down for purposes 
of further calculations. 

To complete the model, a number of emitted photons per 
scintillation pulse must be considered. For this a simple but 
effective and experimentally validated Poisson statistics is em- 
ployed, defining the probability PNin) for the emission of n 
photons; 

-N 

PNin)= — (2.11) 

n'. 

parameterized only by their mean number per scintillation 
pulse. With this final distribution included, a probability density 
Pi(tp, Ie) for first arriving photon being assigned and tp may 
be obtained^: 

oo 

Piitp, Ie) = PN{n)fn{tp, tE\ I) (2.12) 

n=l 

Entering (2.10) and (2.11) into (2.12), while shifting a summa- 
tion index a step backwards; 

oo 

p,(tp, tE) = Ne-''E{tE)A{tp)Si{tp) y -{N[l- I(f£)A,])" 

(2.13) 

an exponential expansion may be recognized; 

Piitp, tE) = A^£(f£)A(fp)5,(fp)e-^'('^'^' (2.14) 

Furthermore, utilizing relation (2.2) yields a probabiUty distri- 
bution DiitA) for the first photon arrival moment f^; 

DiitA) = ^ pi(tp\ tE^tA- tp)dtp (2.15) 

'„„„(') 

Finally, considering a probability to measure the first photon 
arrival time at one end of the fiber at a distance / from scin- 
tillation pulse origin, with arrival time + St at the other end 
distanced L - I, a coincidental timing distribution Ri{6t) for a 
fiber of the length L is found; 

oo 

j D,(tA)DL-i(tA+mtA (2.16) 



From a symmetry in respect to the middle of the fiber, it is ob- 
vious that the following is true"*; 



Rm = RL-i(-6t) 



(2.17) 



The formal grounds for this step are addressed in Appendix B. 
^This is actually on 
dentical for both fib 
discussed in Chapter 4. 



^This is actually only true if the spatial distribution S of emitted photons 
is identical for both fiber ends, i.e. when S{8) = Sin - 6), with S(6) being 



Let us consider the following cases; 
1 .) a photon is emitted at the moment t^ slightly prior to the first 
one reaching the fiber end (f^ < fg), but with a propagation 
time t„ so long that its arrival follows or would follow that 
of the first incoming; 



tE + t„> tE+ tp 



(3.1) 



2.) a photon is emitted slightly later than the first incoming 
(te > Ie)^ but with a propagation time so short that, had 
it not been lost, its arrival at the fiber end would precede the 
first actual one; 

ti- + t„<tE + tp (3.2) 

Evidently, the probability factors and y„, given by (2.6) and 
(2.8), respectively, have to be corrected for such eventualities. 
Therefore, let us denote by Fi the probability for an early- 
emitted photon to actually arrive late at the fiber end. Incor- 
porating the condition (3.1) into integrals' limits and utilizing 
the relation (2.2), the defining expression for Fi becomes; 

ri(tp,tA;l)= J E{t^)dt^ J S i(t„)A(t„)dt„ 

J S i(t„)A(t„)dt„ J E(t^)dt^ 



lA-ts 



(3.3) 

Following the analogue procedure with the probability F2 for a 
late-emitted photon to arrive prior to the defining moment f^, 
the following remains; 

tA-t,„in{l) tA-t^ 

T2{tp,tA\l)= J E(t^)dt^ J S i{t„)A(t„)dt„ 



hniM 



IP 'A -In 

j S ,(t„)Ait„)dt„ j E(t,)dt, 



r,„i„(l) 



tA-tp 



(3.4) 



The equality of two sides in both (3.3) and (3.4) is achieved by 
interchanging the integration order. Depending on the specific 
form of selected distributions, one of the sides may be greatly 
preferred for the actual numerical calculations. 

Since the photons regulated by Fj are, in fact, allowed to 
reach the fiber end, Fi is to be added to the factor a,„: 



am^\ A [I(f£)(l - A;) + Tiitp, tA, I)] 
\m - 1 



m-\ 



(3.5) 



while F2 is to be subtracted from y,„, not allowing any photon 
to precede the first incoming one; 



r,„ = [l-l{tE)-T2{tp,tA\W 



(3.6) 



3 



Following the steps leading from (2.1) to (2.14), it is invari- 
ably found that the corrected form for the probability density 
pi{tp, Ie) equals: 



with a joined corrective probability Viitp, Ia) defined as: 

T,(tp,tA) = r2(tp,tA-J)-rdtp,tA-j) 

which, explicitly written, is of form: 

iA-t„„M 



(3.7) 



(3.8) 



r;(ff,fA)= J EiQdt, j S iit„)Ait„)dt„ 

tA-t,„„Al) t,„M 

AD 



J S i(t„)A(t„)dt„ J E(Qdt^ 



'„uM 



tA-tp 



(3.9) 



An example of rz,/2(fp, Ia) - as shown in Fig. 2 - was calcu- 
lated for a specific selection of parameters and distributions in- 
troduced in Chapter 4. When applied to (2.16), the corrective 
effect from (3.7) would be suppressed by orders of magnitude 
within later presented Figures. Therefore, within Chapter 4 and 
Chapter 6 it was indeed more than sufficient to utilize (2.14) for 
numerical calculations. 
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Figure 2: An example of corrective probability Vi(tp, Ia) for / = L/2. The prop- 
agation time tp is limited by minimal and maximal values ?„„„(/) and f„M.v(0. 
respectively. 



quate form was adopted from [2, 8]: 

1 / fi-'E/T, _ p-tElTl 

EitE) 



r+ 1 



+ — e 

T"3 



-Ie/T3 



@(tE) (4.1) 



with a step function^ ©(?£) regulating causality. A decay time 
Ti - 2.7 ns of the main fast component is reported in [7], while 
the values of additional parameters - T2 - 0.9 ns, T3 - 14.2 ns 
and r = 0.27 - were selected from [2, 8]. 

Path length dispersion S (tp) is deduced from the spatial dis- 
tribution of emitted photons. Assuming isotropic scintillations, 
a normalized probability density S (0) for photon being emitted 
under an angle 6 relative to the fiber axis (Fig. 1) is given by: 



dS(0) = -id(cos6') 



(4.2) 



For the single-cladding fibers, light can be roughly considered 
to propagate only through the scintillating core, since that pass- 
ing through the outermost cladding material is rapidly being 
lost due to the outer surface imperfections and detriments [9]. 
Therefore, eliminating the contribution of the light refracted at 
the core-cladding interface, a photon propagation time tp(6; Z) 
is given by the corresponding path length inside the fiber: 



tp(e;l) 



n I 1 
c COS 9 



(4.3) 



with scintillating core refractive index and the speed of 

light c in vacuum. Utilizing a relation (4.3), (4.2) is translated 
into: 

c \A. ncorel®[tp - tmin(l)\®[tmaxQ) - tp] 

Si{tp)dtp = — ^ dtp (4.4) 

tp 

where f„„„(0 is a minimal time required for the light propaga- 
tion (61 = 0): 

tmM = — (4.5) 

c 

A particular value of ticore - 1-60 is obtained from [7]. 
With a cladding refractive index of tidad - 1.49, tmaxil) is 
a maximal propagation time restricted by the critical angle 
for a total internal reflection off the core-cladding interface 

(cos Oi-i-itical — ndad/f^core)' 



x(/) = 



ndadC 



(4.6) 



4. Model demonstration 

To demonstrate the results of a proposed model, probabil- 
ity distributions with all the required parameters have to be de- 
fined. As the fiber prototype - L - 1.5 m in length - a model 
BCF-10 from Saint-Gobain [7] was selected. Due to the well 
established fact that the simple exponential description of the 
emission probability distribution £(?£) would be an oversimpli- 
fication for the fast emitting plastic scintillators, a more ade- 



Finally, for an attenuation factor a simple exponential depen- 
dency on the path length A(6) - g-'/C'i'^ose) adopted, with a 
single value A = 2.2 m for the attenuation length available [7]. 
Utilizing (4.3), A(ff) translates into: 

A(tp) = e-^'''/(«— (4.7) 



c I VJ 11 ( U 

Throughout the paper, the convention: ©(/) = < . will be used. 
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Figure 3 : Rescaled coincidental timing distributions for 1 = L/2 and tlie average 
ofN= 10' (full line), 10^ (long-dashed line), 10^ (short-dashed line) photons. 
The improvement in resolution for an increased N is evident. 




200 400 600 800 1000 1200 1400 1600 1800 2000 
N 

Figure 4: FWHM dependency on the mean N of emitted photons. The isolated 
frame shows the plot with logarithmic scale for both axes, covering the wider 
range for N. Note that FWHM(A' = lO'') > 0.3 ns. 
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Figure 5: FWHM value dependent on the position of a scintillating pulses' 
point of origin. Full line (left scale) shows the case for N = 4, while dashed one 
(right scale) for A? = lO"*. 



Fig. 3 shows the rescaled coincidental timing distributions 
for I = L/2 and the mean number of N = 10', 10^, 10^ photons, 
emphasizing the decrease in resolution defining FWHM value. 
A fully expected improvement in resolution for the increased 
mean number of emitted photons is explicitly presented in 
Fig. 4. Though, purely mathematically, it is also to be expected 
that FWHM drops to with approaching infinity, even for 
the excessive number of 10^ photons, FWHM amounts to more 
than 0.3 ns. Additionally, for = 4 and = 10"* Fig. 5 shows 
the influence of the scintillating pulses' position along the fiber 
upon the attainable FWHM values. 

In Fig. 6 the coincidental timing distributions are given for 
9 equidistant pulses' origins and separate cases of A^ = 4 and 
A^ = 10"* photons. 



5. Photon count 

Most relevant for the experimental purposes, physical quan- 
tities such as the photon detection efficiency and the mean num- 
ber of detected photons deserve a special consideration. For a 
proposed model (2.16), a photon detection efficiency ^nH) is 
contained within the area below the coincidental timing distri- 
bution curve: 

oo 

*P^(0 - j Rmd(6t) (5.1) 

— oo 

However, due to the complexity of (2.16), little could be 
achieved with the approach from (5.1). Fortunately, there is an 
alternate procedure yielding the exact solution for ^nH)- Let 
us, therefore, consider the probability ?/„(/) for at least one pho- 
ton of total n emitted to reach every end of the fiber; 

«-I n-m / w \ 

- z § pt r>' ^-'^^ - - ^-'^"""'^ ^'-'^ 

Evidently, it is given by all the available combinations of m pho- 
tons reaching one end, k photons the other, with the remainder 
of n-m- ^ being lost. It may be readily shown that ri„(l) equals: 

77„(0 =!-(!- A,)" - (1 - Ai_,)« -f (1 - A, - A^,)« (5.3) 

The transparency of previous result may be appreciated, as it 
excludes the probability for all photons not to arrive at one of 
the ends, while correcting for a double intake of none of them 
reaching either end. This exact form enables a straightforward 
calculation of *I'a'(0^ when the probability for a particular num- 
ber of emitted photons is considered: 

oo 

It is to be noted that (5.4) requires at least two photons, one 
for each fiber end. Prior to disclosing the solution for (5.4), we 
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Figure 6: Coincidental timing distributions for N = 4 (left) and N = Iff* (right). For a given N, 9 equidistant positions of scintillating pulses along the fiber length 
(from right to left) are given by: //L = 0.1i (;' = 1,...,9). 



propose an alternate starting point for calculation of *Pa((Z): 

oo r mix 

11=2 m=l L *=1 ^ ' 

X Pdn - m) r ^ "jAt,(l - Al-/)"-"' 

(5.5) 

based upon the approach discussed in Appendix B and 
considering the arrival of photons at each fiber end separately. 

Both (5.4) and (5.5) consistently yield: 

'i'M = 1 - e-'^'^' - e-^^'--' + e-'^('^'+'^"> (5.6) 



Fig. 7 plots the solution (5.6) for N = 10, 50, 100, 500 pho- 
tons, emphasizing the placement within absolute probability 
scale [0,1]. Fig. 8 shows ^'n(L/2) dependent on average num- 
ber of photons emitted. There is an additional observation to be 
made, based on the fact that for a sufficiently large the detec- 
tion efficiency saturates around 1. It is implied that a limiting 
case of the coincidental timing distribution corresponds with a 
representation of a delta-distribution: 



erage Mn(1): 



Um Ri{6t; N) ^ 6 [6t + t,„in(l) - t„,JL - I)] 



(5.7) 



encompassing, in fact, an infinite set of representations, one for 
every single arbitrary selection of distributions and parameters 
involved. 

Besides the efficiency for detecting any photons at both fiber 
ends, their mean number may also be deduced. For measure- 
ments involving the coincidental timing resolution, a coinci- 
dence of photon arrivals at both ends is a prerequisite that de- 
fines the starting point for the calculation of a coincidental av- 



n-m I \ 



-m-k 



(5.8) 



Due to the coincidence requirement, normalization for the event 
detecting probability less than 1 must be performed, yielding a 
simple solution: 

The fact that (5.9) is devoid of a term Az,_/ is consistent with 
an independence of the fiber ends, formally supported by (B.2) 
( Appendix B). 




Figure 7: Photon detection efficiencies for N = W (short-dashed line), 
50 (medium-dashed line), 100 (long-dashed line), 500 (full line) photons. The 
placement within absolute [0,1] probability range is emphasized. 
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6. Triggering feature 

Photon resolving devices offer an experimental possibility 
for testing the various aspects of the single-photon and/or 
multi -photon statistics. In [3] it is experimentally demonstrated 
that by eliminating the events with low number of detected 
photons (triggering) at either fiber end, the width of the remain- 
ing coincidental timing distribution significantly decreases, i.e. 
resolution improves. Therefore, our goal is to include such 
triggering option within the theoretical model. 

Let us demand that only the events with more than A detected 
photons are acceptable from one of the fiber ends. Hence, the 
probability 7r„,A(0 for at least A out of n photons reaching the 
end is required: 



7r„,A(O = 0(«-A)|]J^jAf(l-A,)« 



(6.1) 



Furthermore, let us consider an effect on a factor p,,, from (2.3), 
which labeling will be extended to /J„,m(0 for further use. Since 
it is assumed within p„j„il) that initial m - 1 of n photons were 
lost, while m-th one does reach the fiber end, n - m photons 
remain for the probability n to regulate: 



Pn,m(.0 



n-m,A 

(I) (6.2) 



Finally, a probability /„ from (2.1) must be corrected for allow- 
ing only the acceptance of more than A photons: 

n 

f„,^(l) = 0(n - A - 1) 2 P„,,„a(0 (6.3) 

m=l 

Upon entering (6.3) into (2.12), the following expression re- 



mams: 



-m-k 



oo n-A n-m i \ 

«=A+1 m=l k=A ^ ' 

(6.4) 

with a triggering parameter A appearing inside the limits of all 
the sums. Evidently, for A = (6.4) reduces to (2.14). Since 
trigger values may be separately set for every fiber end, the final 
coincidental timing model incorporates both of them: 



OO 

J Dl(tA\ 



Ai)DL-i(tA+6t;A2)dtA (6.5) 



Fig. 9 shows the effects of triggering for the case of 
N = 4 photons with Aj = A2 =0, 1, 4, supporting the 
resolution improvement observed in [3]. Moreover, Fig. 10 
presents the reduce in FWHM for N - A and few values of 
Ai = A2 = A. It was found that points from Fig. 10 follow the 
trend FWHM oc (A -i- 1)-" <'24 reasonably well. 

The triggering generalization is also easily extended to the 
model for a photon detection efficiency and a mean number of 
detected photons. Therefore, a trigger inclusive form for the 




100 120 140 160 180 200 
N 



Figure 8: A photon detection efficiency for / = L/2 and a varying mean of 
the photons emitted. An isolated frame presents a logarithmic scale for N, 
encompassing a wider range of values while revealing a probability saturation 
for large N. 



photon detection efficiency equals: 

00 h-I-At 



>I'^(;;Ai,A2)= Z ^^^'"^ Z LK(1-^'> 

«=2+Ai+A2 m=l+A] L /:=1+A| ^ ' 

n-m J \ 



i-k 



-m-k 



(6.6) 



while that for the coincidental mean of detected photons: 
M^(/;Ai,A2)= ^ 



^^(;;Ai,A2) 

M-1-A2 r / \ 



=2+Ai+At m=l+A| 



i:=l+Ai 
n-m 



n-m I X 

PAn-m) Y r:>L(i-AL-,)"-'"-^ 



(6.7) 

By setting all the available trigger values within (6.5), (6.6) 
and/or (6.7) to 0, a model without the triggering effects remains. 
Entering A2 = -1 into (6.7) leaves a non-coincidental average 
of photons detected at the single fiber end. For experiments 
considering not only the photon-inclusive events, but also the 
pedestal measurement, Ai = A2 = -1 may be inserted, simpli- 
fying the photon average calculation: 



Ma,(/;Ai =-1,A2 = -1) = A?A, 



(6.8) 



However, in a course of actual measurements, a single photon 
detection efficiency for the photon resolving detectors must be 
considered. Appendix C addresses this subject, essential for 
the correct estimation of the emitted photons average A^. 
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8t [ns] 

Figure 9: Rescaled coincidental timing distributions for / = L/2 and N = 4, 
showing the improvement in resolution due to triggering. Trigger values were 
set to be equal: Ai = A2 = A, with A = (full line), 1 (long-dashed line), 4 
(short-dashed line). 




Figure 10: Triggering induced decrease in FWHM for / = L/2 and N = 4, 
with Ai = A2 = A. An isolated frame presents both scales as logarith- 
mic, emphasizing the reasonable power-law trend, shown by the full line 

FWHMk(A+ 1)-°-^^*. 



7. Experimental verification 

For the initial verification of the proposed model an ex- 
perimental setup from Fig. 11 was employed, comprising a 
standai-d (Non-S) type SCSF-78 fiber from Kuraray [10] with 
1 mm^ Photonique SSPM_0611BlMM_TO18 [11] silicon 
photomultipliers as central components. The fiber itself is 
= 0.83 mm in outer diameter, with two layers of cladding 
and a spectral maximum of emitted light at 450 nm. The length 
of the selected fiber sample was L - 0.52 m. As a source 
of the scintillating pulses, a yS-radioactive '^^Cs sample was 
used. To ensure their protection from the environmental Ught, 
these elements were enclosed within the plastic container. The 
operating voltage for two diodes of different breakthrough 
values - 26.4 V and 27.0 V - was set to 27.6 V and 28.7 
V, respectively. For a signal amplification the Photonique 
ampHfiers AMP_0604 [11] were used. 

During the final stage of a signal transmission, amplified 
signals were fed to the fast 10-bit 4-channel Acqiris digitizer 
DC282 [12], with the maximal sampling rate of 8 GHz for a 
single active channel. With two active channels required for 
signal intake from both SiPM diodes, signals were recorded 
at 4 GHz sampling rate. Wholly stored for further offline 
analysis, signals may be freely and repeatedly accessed by 
any data extrapolation algorithm. At present time we are 
employing a simple baseline corrected calculation of an area 
below the voltage curve, i.e. signal integration, thus imitating 
the working principle of the charge collecting ADC units. 
For determination of the timing properties, specifically signal 
arriving and rising times, a constant fraction discrimination 
principle of classical TDC units was adopted. For this purpose, 
the leading edge of a signal is fitted to the Gaussian form from 
which the arrival time is extrapolated employing the constant 
fraction factor of fcr = 20%. Fig. 12 presents an example of 
the digitizer recording, together with a leading edge fit. 



At all times during measurements a threshold equivalent to 
a trigger value of Ai = 4 for accepting events from one of the 
diodes was set above the range of a dark noise spectrum in order 
to avoid randomly triggered data collection. Subsequently, all 
the signals from the other diode were indiscriminately kept. Re- 
fined by a reduce in otherwise dominant dark noise due to the 
strict coincidence condition, an energy spectrum of a second 
diode was assumed to follow the multi-Poissonian form <i(x) 
proposed in [3]: 



00 00 



o(v) = y y ^-"n^^-/-^(px)- ^irlJiZt^ ^^^^ 



p=() s={) 



^In [0-2 + {p + i)o-2] 



with X - - xq)/G as a number of fired pixels obtained by 
redefining the original energy scale x in respect to a pedestal 
offset xq and pixel-to-pixel gain G. Summation indices p 
and s enumerate primary and secondary fired pixels^ with 
n and S as their mean, respectively. Parameters ctq and cri 
define the pedestal and pixel noise. Fitting the data, values of 
n = (4.613 + 0.005) and I = (0.467 + 0.002) were obtained. 

For a comparison of the experimental and theoretical results, 
an emission probability distribution £(?£) from (4.1) was 
assumed, with the main component ti - 2.8 ns found in [10]. 
The values of r = 0.27 and T3 - 14.2 ns were kept, while 
the last remaining parameter was set to T2 - 1.1 ns. Though 
in Chapter 4 a fiber model BCF-10 was selected due to its 



'The firing of primary pixels is induced by an actual photon absorption, 
while that of secondary ones by a random dark noise occun'ence, pixel-to-pixel 
crosstalk and after-pulsing [13]. Therefore, with Hp as a mean number of the 
dark noise pulses within a signal integration interval, nc as a crosstalk probabil- 
ity and tta as an after-pulsing probability, the total mean (p + s) of fired pixels 
is given by: 

(p + .v) = n(i + s) = (n + no)(i + ncXi + ^a) 
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Figure 1 1 : Schematic block diagram of the experimental setup for measuring 
the coincidental timing distribution. 



Figure 12: An example of a signal recorded by DC282 digitizer, with discrete 
points sampled every 250 ps. A Gaussian fit of a leading edge is shown in full 
line. 



single cladding, validating the application of a model (4.4) 
for the path length dispersion 5;(fp), the same model was 
retained for SCSF-78 fiber - a practice supported by an inner 
cladding thickness of only 3.4% of a scintillating core diameter 
[10]. Therefore, a light propagation speed inside the fiber 
is considered to be determined only by the core refractive 
index of nc„re - 1-59, while «,„ - 1.49 of inner cladding 
remains omitted within the simplified (4.3) model. However, 
a refractive index n„,„ = 1 .42 of outer cladding is essential for 
the definition of a maximal propagation time t^ax from (4.6). 
Finally, due to the rather short fiber sample used (L = 0.52 
m), a single exponential component of A = 4.21 m [10] was 
adopted for an attenuation factor A{tp). 

To obtain a measurement of the coincidental timing distri- 
bution, a non-coUimated '^^Cs source was placed at the center 
of the fiber: / = L/2. A single photon detection efficiency, 
addressed in Appendix C, was estimated to si = £2 = 20% 
for both SiPM diodes [11]. Based on the value IT ^ 4.6 and the 
calculated photon survival probability A/, a mean number of 
photons emitted during a scintillation decay was evaluated to 
N = 460. 

Fig. 13 presents a coincidental timing distribution remained 
upon triggering with Ai = 4 and A2 = 0, 2, 4, 6. While a trigger 
value Ai was set in a course of the data collection, A2 was 
freely manipulated during offline analysis. Fair, if not excellent 
agreement between the experimental and theoretical results 
may be observed, regarding both the distribution width and a 
shape. 

The photon count considerations brought up in Chapter 5 
must be revised if they are to be translated into a pixel count, 
due to the presence of mechanisms inducing secondary pixel 
discharges, such as the dark noise, pixel-to-pixel crosstalk and 
after-pulsing [13]. When considering only the primary dis- 
charges, i.e. those induced by an actual photon absorption, 
the probability lAn(A) for more than A pixels firing may be ex- 



pressed as [6]: 

^n(A) = 1-2] — (7.2) 

p=() P' 

Therefore, a general trigger inclusive model (6.6) may be ap- 
proximated as: 

'P^(/; Ai, A2) ^ ^n,(Ai)^n,_,(A2) (7.3) 

However, a secondary discharge occurrence somewhat modi- 
fies a purely Poissonian statistics, increasing the probability for 
pixel activation, extending (7.2) into: 

p=i) j=o ' ■ 

Fig. 14 plots the experimentally obtained decrease in a number 
of registered events against an approximate expression (7.4), 
together with the result of a secondary-discharge-free model 
(7.2). 

8. Conclusions 

The scintillating fibers provide an efficient experimental 
tool for the detection of charged particles and at least a partial 
reconstruction of their trajectories. Due to the reconstructive 
methods relying upon the time of propagation of the scintil- 
lation light, timing properties of the scintillating fibers have 
become a central issue of their experimental performance. 
Therefore, understanding the sources of their clearly quan- 
tifiable coincidental timing resolution is a requirement of a 
fundamental importance for further experimental applications. 
From a developed model it was found that the rate of a scin- 
tillation decay, the number and spatial dispersion of emitted 
photons are defining factors for such resolution. Additionally, 
it was confirmed that resolution may indeed be improved by a 
triggering procedure, though implying the penalty of a reduced 
photon detection efficiency. Experimental results support 
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the developed model, indicating, however, the necessity for 
a clear identification of a portion of random events within 
performed measurements. This is achieved by considering a 
physical impact of the relevant experimental elements besides 
the scintillating fiber itself, e.g. photon resolving devices. 
Further experimental investigation is planned in order to de- 
termine the practical applicability and limitations of such setup. 
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Appendix A. Equivalency proof 



Given the general emission time probability distribution fit) 
- previously denoted as - a probability distribution e{t„) 
for the emission of the first of total n photons at the moment 
/„ is to be determined. There are two basic approaches to this 
procedure - with or without arranging all the n photons in time. 
A model including the time arrangement is defined by: 

OO DO OO 

E(t„) ^ n\f(t„) J f(t„-i)dt„^i J f(t„-2)dt„-2--- J /(fi)dfi 

tn 'n-l h 

(A.l) 

with the factorial term determining the number of photon per- 
mutations. Term /(f„) is a probability density for the emission 
of first photon at the moment f„, while every subsequent inte- 
gral regulates the probability for a given photon to be emitted 
following the previous one (f, > f,_i for (n + I - /)-th photon). 
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Figure 14: Triggering induced decrease in a number of registered events, with a 
constant value of Ai = 4. Experimental points are shown as circles. Rectangle- 
shaped points are the result of an approximation (7.4), while the triangular ones 
of a model (7.2) devoid of the secondary pixel discharges. Dashed lines serve 
as a guide for an eye. 



On the other hand, a model without arrangement requires: 

OO OO DO 

E(tn) = nf(t„) J f(t„-i)dt„^i J f(t„-2)dt„-2 ■ ■ ■ J f(h)dh 

(A.2) 

where initial n regulates the number of possibilities to select 
only one (fist) photon, with all remaining n - 1 following in 
any order Hence, the difference between (A.l) and (A.2) is 
contained within combinatory factor and lower limits of inte- 
gration. Equivalency between (A.l) and (A.2) would enable us 
to confidently and without reservations use whichever model 
proves to be more convenient for further calculations. 

Thus, we begin the equivalency proof by considering the 
properties of a basic primitive function F{t) for the emission 
probability distribution /(f). The following holds: 



F(oo) = 



(A.3) 



This statement is supported by a physical requirement for 
the emission process to have a well defined mean lifetime 
(t) = tf{t)dt. For this to be true, it is obvious that /(f) is 
required to asymptotically decrease strictly faster than 1/f^, di- 
rectly implying that the primitive function F{t) decreases faster 
than 1//, validating the claim (A.3). Moreover, from the proba- 
bility normalization: 



= J f(t)dt^F(oc)-F(0) 





directly follows: 



F(0) = -1 



(A.4) 



(A.5) 



Having obtained the result (A.3), it is easily shown that all the 
identical integrals in (A.2) give rise to the following solution for 



e(f„): 

s(tn) = n/(f«) 



OO 

/ 



/(f)df 



-(-l)"-'«/(f„)^'""'(f„) (A.6) 



In order to prove that a model (A.l) with photon arrangement 
yields the identical solution, a method of mathematical induc- 
tion is applied. For this purpose we define the term I„(t„) as: 

OO OO OO 

In(tn)= J /(f«-l)df„-l J f(t„-2)dt„-2--- J /(fl)dfl (A.7) 



and state that the following is valid: 

In(t„) 



(n-1)! 



F"-\t„) 



(A.8) 



Previous expression will serve as an induction hypothesis, sup- 
ported by the basis: 

/i(fi) = 1 (A.9) 

Within the inductive step it is to be showed that the hypothesis 
holds for/„+i(f„+i): 



4+i(f«+i) 



OO 

J f(t„)In(t„)dt„ 



(A. 10) 



Upon entering (A.8) into (A. 10) and partially integrating: 

OO 



4+l(fn+l) 



(n-1) 



J f(t„)F"-\tn)dt„ 



u^F"-\t„) ^ du^(n-l)F"-Ht„)f(t„)dt„ 
dv = /(f„)df„ ^ V = F(tn) 



i-iy 



n-1 



(n - 1)! 



F"(tn)\Z,-(n-l) j f{t„)F"-\t„)dt„ 
',,+1 

(A. 11) 



the appearance of the same integral may be noticed on both 
sides of the equation. Rearranging the terms leaves: 

4+i(f„+i) = ^^^'"(Wi) (A. 12) 

n! 

concluding the inductive proof. 

Finally, the expression (A.l) is identified as a: 

e(f„) = n!/(f„)/„(f„) (A. 13) 

so that introducing (A.8): 

s{t„) = (-l)"-'«/(f„)F"-'(f„) (A.14) 
we may witness a result identical to (A.6). 
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Additionally, with only a step further it is shown that the 
probability distribution s{t„) is, indeed, normalized, since: 

oo 

J £(r„)df„ = «!/„+i(0) = (-irF"(0) = 1 (A.15) 



Here the identity (A. 5) was used. 

Appendix B. Model justification 

Backtracking the origin of the expression (2.16) to (2.12), a 
particular explicit assumption of the proposed model becomes 
evident - the numbers of photons (denoted by n) emitted to- 
wards the two ends of the fiber are considered independent, with 
a mean number of total emitted apparently being identified 
with a mean for every fiber end. Arguably, the more transparent 
approach would be to simultaneously consider both fiber ends 
during the emission of total n photons, with as their genuine 
mean, while distributing them between those ends: 

oo n-l 

= Z Z PN(m)Ml)PN{n - m)fn-,n{L - (B.l) 

It is to be noted that the probability density P/ for arrival of 
photons at both fiber ends requires a minimum of two photons 
to be emitted, while leaving at least one of them for every end. 
Upon entering (2.10) and (2.11) into (B.l), it may be readily 
shown that the following holds: 

P/ = PiPL-i (B.2) 

which is the formal ground for the separate and independent 
treatment of fiber ends. 

Appendix C. Implementing advanced physical models 

Considerations contained within (4.3) and (4.6) are basically 
equivalent to the meridional approximation, which does not 
provide a complete description for the geometry of photon 
propagation [9, 14]. However, developed for a fully general 
selection of relevant distributions and parameters, a model 
comprised within (2.16) is readily available for incorporating 
any improvements and/or additions to the physical description 
from Chapter 4. Therefore, it is our goal to propose the manner 
in which to include a wide variety of advanced physical con- 
siderations. A more detailed overview of these contributions 
may be the subject of future work. 

To improve upon the meridional approximation itself, a prop- 
agation time tp in (4.3) should be adequately redefined in terms 
of additional parameters, e.g. a skew angle, alongside primary 
9. Furthermore, such redefinition enables the inclusion of the 
multiple claddings via separation of the regimes for photons 
trapped inside a scintillating core from those passing through 



one or more claddings. Therefore, given a set {0**') of the rele- 
vant parameters, a propagation time should be defined as: 

tp = in n ® {^f^ - ^'^) ® - ^s) 

I k 

(C.l) 

with 0^*^' as a boundary of an /-th separate parameter space 
domain for the parameter 6*'*'. Evidently, within (4.3) it is 
assumed: 0^'^ = 9, with 0^ ' = and - 9cnticai- In order not 
to omit a potential difference in the attenuation mode within 
separate fiber layers, an attenuation factor A{tp) should be 
redefined in parallel with tp, in a manner identical to (C. 1). 

The following consideration brought up in [9] is a Fresnel 
reflection above the critical angle for a total reflection. This 
effect is easily included by expanding the path length disper- 
sion from (4.4) beyond previously imposed sharp cut, while 
setting tmax — > oo. At this point it is to be noted that path length 
dispersion S i{tp) inherits the normalization from its spatial 
origin, i.e. (4.2), making the integral ^ S i(tp)dtp not equal to 
unity, which is crucial for normalizing (2.16). 

The finite transmittance, i.e. reflectance q less than 1 even in 
a course of the otherwise total reflection, is to be implemented 
through the attenuation factor: 

A(tp;q)^q''^''-'>A(tp;q^l) (C.2) 

with K(tp) as the number of reflections off" the inner fiber in- 
terfaces [9, 14]. Reflections off the fiber ends present an addi- 
tional matter that may be considered when describing the statis- 
tics of photon propagation. For this purpose an extended path 
length dispersion S (tp; I) may be constructed using the original 
reflections-free form 5;(fp): 

OO 

S{tp; = (1 - 2) 2 e^*: YS,^2,i^{tp) + QS 2(k+l)L-l(tp)] (C.3) 
k=0 

to enumerate the events with an even and odd number of reflec- 
tions off the fiber ends prior to the final photon transmission. 
The reflectance Q (Q + q) determines the probability for a 
given number of reflections as well as the final transmission 
probability (\-Q). Evidently, for 2 = an extended dispersion 
S(tp; I) reduces to S i(tp). 

Many, if not all physical properties display a strong depen- 
dency on the emitted light wavelength A - from the photon 
emission itself, to refractive indices of the fiber, attenuation 
lengths, etc. Given the normalized wavelength emission spec- 
trum p{A), this dependency is to be included by averaging the 
probability density p(tp, f^) from within (2.15): 

p(tp,tE)^ j p(tp,tEU)p(m^ (C.4) 

over the relevant wavelength interval extending from /!,„,„ to 
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included in a simple manner: 



A(tp; Si) = SiA(tp; ei = 1) 



(C.8) 



During the charged particle's transition through the fiber ma- 
terial, an entire array of scintillating molecules is excited along 
its path, which generally may not be assigned a unique value for 
the distance / from one of the fiber ends. Thought the basic prin- 
ciple for including photon origin variations remains the same as 
for wavelength variations, it requires somewhat more complex 
approach. Therefore, let us assume that for every single particle 
passing through the fiber, their trajectory is parameterized by X 
- a single parameter or, symboUcally, a set of parameters. For a 
collimated particle beam this parameterization may be omitted. 
However, its significance becomes apparent, for example, in a 
case of a non-coUimated isotropically emitting particle source. 
When including this improvement in respect to the point-like 
description of scintillating pulses, the first step is to translate a 
probability distribution Z)K^a) from (2.15) into a form !Z)(?a; X): 

'D{tA\X)= j Dnx){tA-tiiX))Al{X) (C.5) 

taking into account that / varies from l„,i„(X) to l,„ax{X)- There 
is an additional subtle consideration to be noted - the inclu- 
sion of the moment tux) of a particle reaching the distance 1{X) 
from one of the fiber ends. Its definition may be as simple as: 
hix) = [liX) - Iq{X)\Ivi with Iq{X) as a value for 1{X) at the mo- 
ment of particle entering the fiber and v; as its speed along the 
fiber axis. The following step: 



Furthermore, in a course of the signal timing measurements a 
spread critA) due to any kind of fluctuations within a light de- 
tecting system - electronic noise, fluctuations during the elec- 
tric charge collection, variations in a rise time of a signal, etc. 
- may be taken into account. Consequently, it is manifested 
through the modified arrival times distribution ^;(^a): 



SliitA) = (Di * o-)itA) = Di{i)(r{tA - t)dt 



(C.9) 



gained by means of a functional convolution. 
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!USt;X) 



OO 

-I 



D{tA\X)'D{tA + 6t\X)d.tA 



(C.6) 



brings nothing new to the practice from (2. 16), save for inher- 
iting the X-dependency from 'D{tA;X). However, it constitutes 
only an intermediate step preceeding the final calculation of a 
coincidental timing distribution R(i){6t): 



R(i)m 



ia6t;X)p{X)dX 



(C.7) 



which is - upon averaging procedure considering a distribution 
p{X) - parameterized by a mean value (/). At this point an 
important note has to be made regarding an integration being 
applied to D,{tA) in (C.5), while in (C.7) to 3?,(5f;X). The 
practice from (C.5) allows the photons reaching separate fiber 
ends to be assigned different origins 1{X), with a practice from 
(C.7) compelling all the photons from a single scintillation 
pulse to be assigned a single particle trajectory. 



In the end, even the properties of the light detection devices 
may be implemented through a proposed model. A single pho- 
ton detection efliciency e, for a detector at j-th fiber end is to be 
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